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ABSTRACT 

We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The 
model is based on the concept of a relativistic jet interacting with compact gas condensations produced when 
a star (a red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of 
the evolution of the envelope lost by the star within the jet, and calculations of the related high-energy radiation. 
The model readily explains the day-long that varies on timescales of hours, GeV gamma-ray flare from 3C454.3, 
observed during 2010 November on top of a plateau lasting weeks. In the proposed scenario, the plateau state is 
caused by a strong wind generated by the heating of the stellar atmosphere due to nonthermal particles accelerated 
at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant 
after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key 
model parameters of the source, including the mass of the central black hole: Mbh — 10 9 M Q , the total jet power: 
Lj ~ 10 48 erg s -1 , and the Doppler factor of the gamma-ray emitting clouds: 8 — 20. Whereas we do not specify 
the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in 
the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared 
to other radiation channels of directlyaccelerated electrons. An injected proton distribution oc E~ 1 or harder below 
the relevant energies would be favored to alleviate the tight energetic constraints and to avoid the violation of the 
observational low-energy constraints. 
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1. INTRODUCTION 

3C454.3 is a powerful flat-spectrum radio quasar located at 

a redshift c rs = 0.859. This source is very bright in the GeV 

energy range; during strong flares, its apparent (isotropic) lumi- 
nosity can reach Ly > 10 5 ° erg s 1 (e.g., Striani et al. 2010; 
Ackermann et al. 2010; Vercellone et al. 201 1 ; Abdo et al. 201 1). 
The mass of the central black hole (BH) in 3C454.3 is estimated 

to be in the range M BH ~ (0.5-4) x 10 9 M G (Gu et al. 2001; 
Bonnoli et al. 2011). This mass implies an Eddington lumi- 
nosity Lem ^ (0.6-5) x 10 47 erg s _1 , which is several orders 

of magnitude below L y . Although the large gap between L Hdd 

and L y is naturally explained by relativistic Doppler boosting, 
estimates of the jet power during these flares appear, in any 
realistic scenario, close to or even larger than the Eddington 
luminosity (Bonnoli et al. 201 1). Even though they are quite ex- 
treme, “super-Eddington” jets nevertheless cannot be excluded 
for accreting BHs with very high accretion rates but low ra- 
diation efficiencies. Although 3C454.3 is an exceptional case 
(see, for example, the discussion in Bonnoli et al. 2011), data 
from other objects provide additional evidence in favor of jets 
with super-Eddington mechanical powers (Lopez-Corredoira & 
Perucho 2012). 

The GeV emission from 3C454.3 is highly erratic, with 
variability timescales as short as 3 hr, as reported, in particular, 
for the 2009 December flare (Ackermann et al. 2010). The most 

spectacular flare in terms of both variability and gamma-ray 


luminosity was observed in 2010 November by the AGILE 
and Fermi Large Area Telescope (Fermi /L AT: Vercellone et al. 
201 1 ; Abdo et al. 2011) telescopes. During this high state, with 
the most active phase lasting for five days, the apparent GeV 
luminosity achieved was L 7 « 2 x 10 5 * * * ° erg s' 1 . Around the 
flare maximum, the rising time was t r ~ 4.5 hr and the decay 
time was tf 15 hr. The detection of photons with energies up to 
~30 GeV, the short variability, and the contemporaneous X-ray 
flux constrain the Doppler boosting of the emitter to <5 m j„ >16 
to avoid severe internal yy -absorption in the X-ray radiation 
field (Abdo et al. 2011). 

A remarkable feature of the gamma-ray emission from 
3C454.3 is the so-called “plateau phase” revealed during the 
bright flare in 2010. This plateau phase is characterized by a 
long-term brightening of the source, a few weeks before the 
appearance of the main flare. Such plateau states have been 
observed by Fermi /LAI for three flares (e.g., Ackermann et al. 
2010; Abdo et al. 201 1), with the plateau emission being about 
an order of magnitude fainter than that of the main flare. 

Remarkably, the rapid gamma-ray variability of 3C454.3 is 
accompanied by activity at lower energies. Simultaneous multi- 
wavelength observations of the source during flares have re- 
vealed a strong correlation between optical and X-ray data. This 
correlation has been interpreted as evidence that the gamma-ray 
source is located upstream from the core of the 43 GHz radio 
source, which is at a distance z < few pc from the central BH 
(see, e.g., Jorstad et al. 2010, 2012; Wehrle et al. 2012). 
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Over recent years, several works have attempted to explain 
the flaring gamma-ray activity of 3C454.3 within the framework 
of the standard synchrotron self-Compton (SSC) or external 
inverse-Compton (EIC) models (Katarzynski & Ghisellini 2007; 
Ghisellini et al. 2007; Sikora et al. 2008; Bonnoli et al. 201 1). 
In the SSC scenario, it is possible to reproduce the spectral 
energy distribution (SED) from optical wavelengths to gamma 
rays. In these models, most of the jet power is (unavoidably) 
carried by protons, and only a small fraction is contained 
in relativistic electrons and the magnetic field. The required 
proton-to-Poynting flux ratio, L p /L B ~ 100, is quite large. 
Such a configuration would be hard to reconcile, at least in 
the gamma-ray emitting region close to the central BH, with 
an undisturbed jet that is launched by the Blandford-Znajek 
(Blandford & Znajek 1977) type (BZ) process, in which the 
luminosity of the jet is dominated by Poynting flux and the jet 
consists of e ± -pairs (see also Ruffini & Wilson 1975; Lovelace 
1976). In this regard, we should mention that recent relativistic 
magnetohydrodynamical (MHD) simulations of jet acceleration 
yield much less efficient conversions of magnetic energy into 
bulk motion kinetic energy; these calculations (Komissarov et al. 
2007, 2009; Tchekhovskoy et al. 2010) predict a quite modest 
ratio (L p + L e ±)/L B < 4. 

As stars and clumpy matter are expected to be present in the 
jet surroundings, they could be behind the powerful gamma-ray 
events in active galactic nuclei (AGNs; see, e.g., a discussion 
in Bosch-Ramon et al. 2012). In particular, a red giant (RG) 
can enter into the jet, lose its external layers, and in this way 
generate a strong perturbation inside the jet. This perturbation 
can accelerate particles and produce high-energy radiation. 

The jet-RG interaction (IRGI) scenario has been invoked 
to explain the day-scale flares in the nearby non-blazar-type 
AGN M87 (Barkov et al. 2010; Blandford & Konigl 1979). The 
JRGI has also been applied to the TeV blazar PKS 2155—304 
(Barkov et al. 2012a) to demonstrate that the jet-driven acceler- 
ation of debris from the RG atmosphere can explain the ultra- 
fast variability of very high-energy (VHE) gamma-ray emis- 
sion on timescales as short as r ~ 200 s. A distinct feature of 
the JRGI scenario is the high magnetization (L B /L p e 1) 
of the relativistic flows located at sub-parsec distances, where 
the gamma-ray production supposedly takes place. Although 
the strong magnetic field, B ^ 10 G, dramatically reduces the 
efficiency of the inverse Compton scattering of electrons, it 
opens an alternative channel of gamma-ray production through 
synchrotron radiation of protons (Aharonian 2000; Miicke & 
Protheroe 2001). The latter can be effectively realized only 
in the case of acceleration of protons to the highest possi- 
ble energies, up to 10 20 eV. Thus, the second (somewhat “hid- 
den”) requirement of this model is a very effective accelera- 
tion of protons with a rate close to the theoretical limit dic- 
tated by classical electrodynamics (Aharonian et al. 2002). Al- 
though somewhat speculative, this condition does agree with 
rather model-independent (derived from first principles) argu- 
ments that the relativistic jets in AGNs present the best can- 
didate sites of production of the highest energy cosmic rays 
(Aharonian et al. 2002; Lemoine & Waxman 2009). The large 
magnetic field and acceleration of protons to the highest possi- 
ble energies, coupled with strong Doppler boosting in relativistic 
outflows, not only provide an extension of the gamma-ray spec- 
tra up to TeV energies, but also produce variability timescales 
as short as 1 hr (Aharonian 2000). 

It is interesting to also note that inverse Compton models can 
be accommodated, at least in principle, in the JRGI scenario. 
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Moreover, unlike SSC models of powerful blazars, in which 
the requirement of a very weak magnetic field (implying 
a deviation from the equipartition condition by orders of 
magnitude) generally is not addressed and explained, the JRGI 
scenario can offer a natural way for leptonic models to be 
effective assuming that the gamma-ray emission is produced 
through inverse Compton scattering in shocked clouds originally 
weakly magnetized (see Barkov et al. 2012b). 

In this work, we show that the JRGI scenario is a viable 
mechanism for the explanation of the flares seen in 3C454.3. 
We also argue that within this model, the plateau state can form 
due to the interaction of the jet with a stellar wind excited by 
nonthermal (accelerated) particles that penetrate into the RG 
atmosphere. Given its extreme nature and importance, we will 
use the very powerful GeV flare of 3C454.3 that occurred in 
2010 November as a template for our interpretations. We will 
show that this active period of 3C454.3, consisting of day-long 
flares with variabilities as short as a few hours on top of a plateau 
state lasting several weeks, can be explained by the interaction 
of a RG with the jet relatively close to the central BH. 

2. THE STAR-JET INTERACTION SCENARIO 

The possibility that some exceptional flares in AGNs may 
originate from star-jet interactions has been proposed and 
discussed in our recent papers (Barkov et al. 2010, 2012a, 
2012b). In this work, we explore the question whether the 
extremely bright and short gamma-ray flares detected by Fermi/ 
LAT from 3C454.3 in 2010 November can be explained by 
this model. It is important to note that these flares provide us 
with information of exceptional quality concerning both the 
temporal behavior of the source and its energy spectrum. In 
this regard, we should note that the operation of Fermi /LAT in 
scanning mode allows continuous monitoring of the source, so 
its temporal behavior, even during the pre- and post-flare epochs, 
can be studied without interruption. Before considering the 
case of 3C454.3, we discuss in this section the different 
implications of the JRGI scenario regarding the general temporal 
structure of the active phase in a blazar-type AGN. 

When a star enters an AGN jet, some stellar material is ex- 
pected to be released inside the jet flow, eventually forming 
a population of blobs that will be accelerated by the jet ram 
pressure. During this process, depending on the local magnetic 
fields, different types of shocks and plasma waves can be pro- 
duced and propagate through both the jet material and the blobs. 
These waves may accelerate particles to relativistic energies. 
The interaction of these particles with the ambient radiation and 
magnetic fields would result in broadband nonthermal emission. 
Generally, the properties of this emission depend on the specifics 
of the acceleration and radiation mechanisms and the details of 
the target fields. However, in the case of extremely effective 
particle acceleration and radiation, i.e., when the corresponding 
acceleration and radiative cooling times are shorter than other 
timescales characterizing the dynamics of the system (a manda- 
tory condition given the enormous luminosity of the emission), 
the lightcurve of the emission will be essentially determined by 
the jet/blob interaction model. 

In the fast cooling regime, the proper intensity of the nonther- 
mal emission, i.e., the intensity in the blob comoving reference 
frame, is proportional to the energy released at the jet-blob 
interface. This energy release can be characterized by a sim- 
ple dynamical model, which describes the acceleration of the 
blob by the jet ram pressure. In this model, there are just a few 
relevant parameters that describe the basic properties of the jet 
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and the blob: the jet ram pressure (Pj) and bulk Lorentz factor 
(Fj) and the blob mass (M b ) and radius (r b ; or, equivalently, its 
cross-section: S b = n for details, see Barkov et al. 2012a). 
The time dependence of the intensity of the jet/blob interaction 
corrected for Doppler boosting can be treated as a first-order 
approximation for the radiation lightcurve. 

There are several important stages in the JRGI scenario that 
may have an important impact on the observed lightcurve of 
the radiation: (1) first, the removal of the star’s external layers; 
(2) the formation of a cloud from the removed stellar material, 
which under the jet impact suffers heating and expansion, and 
eventually fragments into a set of smaller blobs 8 (see, e.g., 
Barkov et al. 2012a; Bosch-Ramon et al. 2012; Imshennik 
1972); (3) the acceleration of the blobs in the direction of the jet 
flow to a Lorentz factor equal to or smaller than the jet Lorentz 
factor. The intensity of the nonthermal processes associated with 
the blob motion has a strong time dependence. At the start of the 
acceleration of the blob, its Lorentz factor is modest, yielding 
a small (but growing with time) boosting of the flux. At later 
stages, when the blob Lorentz factor approaches the Lorenz 
factor of the jet, the radiation intensity decreases because of 
the weakening of the jet-blob interaction. According to the 
results obtained in Barkov et al. (2012a) in the fast cooling 
limit, the maximum of the apparent luminosity, accounting for 
the Doppler boosting of the emission, occurs when the blob bulk 
Lorentz factor reaches T b ~ 0.8Fj. 

During the penetration of the star into the jet, strong shocks 
are induced in the stellar atmosphere and a significant part of the 
envelope can be removed. Although the details of this process 
might be very complicated and depend on the radial density 
profile of the star and the structure of the jet outer boundary, 
here, for simplicity, we assume an instantaneous penetration. 
To a certain extent, this simplification is justified by the stellar 
proper motion, which is expected to be faster than the initially 
induced shock waves. Nevertheless, independent of the details 
of the star penetration, one can expect that when the star is fully 
within the jet, a dense cloud will be released. Later, as the cloud 
propagates through the jet, the star may still release an intense 
wind or a sequence of small clouds due to ablation of the stellar 
atmosphere by the jet. A sketch of the JRGI scenario is presented 
in Figure 1 . 

The mass of the cloud AM formed at the initial stage can be 
estimated by comparing the jet ram pressure, 

Lj 
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Figure 1. Sketch of the JRGI scenario, in which a star moving from left to right 
penetrates the jet. The star’s external layers are shocked and carried away, and 
a cometary tail, the origin of the plateau emission, forms. The acceleration and 
expansion of the bigger clouds from the initially blown-up external layers of the 
star would lead to the main flare. 


(A color version of this figure is available in the online journal.) 


the star is 


M = 


cttPqRI 

CM, 


( 3 ) 


A more precise estimate of this rate requires a detailed study 
of the interaction process, which is beyond the scope of this 
paper. Nevertheless, in Section 3.2 we discuss whether a 
sufficiently strong stellar wind can be excited by the penetration 
of nonthermal (accelerated) particles. 

Since the initial size of the expelled cloud should be compa- 
rable to the size of the star, it is possible to estimate the cloud 
expansion time as f exp oc 2 R*/c s , where c s is the sound speed of 

the shocked material: c s ~ [(47rP 8 /3)y g .Po/M c ] 1 ^. The cloud 
expansion time is 


‘exp 


*exp 


M c 

y g R*Pp 


1/2 


( 4 ) 


where y g — 4/3 is the plasma adiabatic coefficient and A exp is 
a constant of about a few (Gregori et al. 2000; Nakamura et al. 
2006; Pittard et al. 2010; Bosch-Ramon et al. 2012). According 
to the relativistic hydrodynamic (RHD) simulation by Bosch- 
Ramon et al. (2012), a value of 1.5 can be adopted for A exp . 

The blob acceleration occurs on a timescale of (see, e.g., 
Barkov et al. 2012a): 


with the gravitational force. Here, cop is the jet cross-section at 
the star crossing height. This gives the following estimate: 


a if D < 1 

al if D > 1. 

c D 


( 5 ) 


AM = 


xPoRj 
GM * ’ 


( 2 ) 


where M* and R, are the RG mass and radius, respectively. 

The absolute upper limit on the continuous mass-loss rate 
can be estimated by comparing the energy flux density of the 
jet qo = LJttcoq and the gravitational force work density for 
removing matter from the stellar surface q hd = pvG M t J R t . = 
MGM*/ [j: Rl). The maximum continuous mass ejection from 


8 Throughout the paper, we will use the subscripts “c,” “b,” and “0” to 
characterize the parameters of the initially formed cloud, of the blob, and of 
the jet at the star crossing height, respectively. 


The D-parameter that will be often used in the paper has a simple 
meaning. It is the dimensionless inverse mass of the blob: 


PpJtr^zp 
4c 2 M b T 8 ' 


(6) 


The above timescale corresponds to the blob acceleration in 
the laboratory reference frame. However, since the blob gets 
accelerated toward the observer, the emission delay, as seen by 
the observer, should be approximately corrected by a factor of 
l/(2To). Thus, the observed peak of the emission should be 
delayed by a time interval of 


A t — ? exp + Lee/ (2T /) . 


( 7 ) 
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There are some uncertainties in this equation. In particular, for 
D < 1, the blobs can travel a distance comparable to z and 
the jet properties at the dominant emission blob location may 
differ significantly from those at the blob formation site. Thus, 
the detailed evolution of the emitter can only be addressed 
properly through numerical modeling. However, despite the 
adopted simple approach, some important conclusions can be 
derived. Namely, given the dependence of the two components 
on the right-hand side of Equation (7) on the stellar material 
mass ( M c and M b , respectively), the radiation associated with 
the heavy cloud expelled first will be delayed with respect to the 
emission produced by lighter clouds formed later. 

The emission produced by lighter clouds allows an estimation 
of the time required for the star to cross the jet. Once the star 
enters into the jet, the process of jet-star interaction should 
proceed steadily, with the production of these lighter clouds 
being roughly constant on average. Thus, the whole duration of 
the light cloud-associated emission, if observed, can be taken 
as a direct measurement of the jet crossing time to 2a>o/ Vjnb, 
where V orb < yJ2GM^/zo is the stellar velocity. This yields 
the following relation: 


to > f (8 ) 

cr g 

where r g is the BH gravitational radius. Adopting the paradigm 
of magnetically accelerated jets (see Appendix A, Equa- 
tions (A2)-(A4)), it is possible to derive a very simple expression 
for this timescale: 

to > 2 3/2 zo/c . (9) 

In this way, the duration of the jet-star interaction is determined 
by the interaction distance from the central BH. 

The physical properties of the emitting blobs determine the 
available energetics and the maximum flux reachable in the 
considered scenario. In the blob comoving frame, the jet-star 
interface energy flux is defined by 


qh [n r 4 ] 4 • 


( 10 ) 


Assuming a fixed efficiency £ in the blob reference frame 
for the transfer of jet power to nonthermal gamma rays (where 
<gC 1), and correcting for Doppler boosting, one can estimate 
the luminosity of a blob: 


which accounts for the total energy transferred by the jet to 
a blob during the acceleration process, MbC 2 ro, and for the 
anisotropy of the emission due to relativistic effects represented 
by the factor Tg (see Appendix B, in particular Equation (B9), 
for details). 


3. THE 2010 NOVEMBER FLARE 
3.1. General Structure of the Active Phase 

The total apparent energy of the GeV gamma-ray radiation 
detected during the flare observed from 3C454.3 in 2010 
November was about E tot ~ L y At/(l + Zis ) ~ 2.3 x 10 55 erg. 
The exceptionally high flux during this period allows the 
derivation of a very detailed lightcurve, as seen from Figure 1 
in Abdo et al. (201 1). The nonthermal activity lasted for ff u n ~ 
80 days. The onset of the activity period was characterized by a 
plateau stage. During the first f p i ~ 13 days, a rather steady flux 
was detected, with an apparent luminosity L p i ~ 10 49 erg s -1 . 
The plateau stage was followed by an exceptionally bright 
flare, the total duration of which was fa ~ 5 days, with a 
rise time of f r ~ 4.5 hr. The maximum flux reached was 
7 x 10~ 5 photons cm 2 s -1 , which corresponds to a luminosity 
of Ly ~ 2 x 10 5 ° erg s _1 . The final stage of the flare phase was 
characterized by variable emission with a flux approximately a 
factor of ~5 weaker than the main flare, but still a factor of ~2 
above the plateau level. 

The observed luminosity of the plateau phase allows us 
to determine a lower limit on the star mass-loss rate by 
differentiating Equation (12): 

M*^ lO^Lpj^r'r^gs- 1 , (13) 

where L p i 49 = L p i/ 1 0 49 erg s' 1 . Comparing this requirement 
to Equation (3), it is easy to see that for typical parameters of RG 
stars, Equation ( 1 3) represents a modest fraction (a few percent) 
of the maximum possible mass-loss rate. In Section 3.2, we will 
return to this issue. 

To evaluate the feasibility of the IRGI scenario for the 
3C454.3 main flare, it is necessary to check whether the flux, the 
total energy release, and the flare delay with respect to the onset 
of the plateau are well described by Equations (7), (11), and (12) 
for a reasonable choice of jet/star properties. For example, 
Equation (11) can be rewritten as 

Po = 8 x 10 3 Ff^LyjoT^ ergcnT 3 , (14) 


Ly =^cF e P 0 rljt,i, (11) 

where the correction function F e depends on time, or, equiva- 
lently, on the blob location in the jet (for a mathematical def- 
inition of this function, see Appendix B, Equation (B8)). We 
note that the structure of the jet, i.e., the dependence of the jet 
Lorentz factor on z, determines the actual dependence of F e on 
Z (also see Appendix B for details). 

The maximum value of F e monotonically depends on the 
/9-parameter, approaching a value of 0.4 if I) > 1 and a value 
of ~0.1 for D = 0.1. This relatively weak dependence allows 
us to derive the maximum blob luminosity. Also, it is possible 
to obtain an estimate of the total energy emitted by a blob or 
an ensemble of sub-blobs as a result of the fragmentation of the 
original cloud (M c = Mb): 


where L y so is the peak luminosity of the flare normalized to 
10 5 ° erg s -1 and .S’1-, 32 = ztr^/ 10 32 cm 2 is the blob cross-section. 

A total energy budget of the flare of ~2 x 10 55 erg is feasible, 
according to Equation (12), if 


M c .mT; 


C.301 0,1.5 


0.04£y, ss 

H Pe 


0.1 

We’ 


(15) 


where M Ci 30 = M c /10 30 g is the mass of the blown-up RG 
envelope (i.e., the initially formed cloud). This requirement 
appears to be very close to that provided by Equation (2), where 
the jet ram pressure is now substituted using Equation (14): 


Mr- 


5 x 10 29 


^*,2^*,0^y.50ro,1.5^b,32 S' 


(16) 


E y ~ 8£F e M b/c c 2 r 3 , 


(12) where = R*/10 2 Rq and M* o = M* / M 0 , respectively. 
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The second term in Equation (7) is expected to be short 
compared to the duration of the plateau phase, even for D ~ 0.1, 
and thus the duration of the initial plateau phase constrains the 
expansion time (see Equation (4)): 

i exp » 5.4 x UfFl^j'r- 

x <M /2 C-^so r 0l i.5<^ S . (17) 


where z is the distance from the central BH. It can be seen from a 
comparison of Equation (20) and Equation (24) that the solution 
found implies a jet with a luminosity exceeding the Eddington 
limit by a factor of > 10. This requirement is also consistent 
with the lower limit on the jet luminosity 

Lj > cS h P 0 = 8 x 10 47 Zq i 7 § -1 ^“ ergs -1 , (25) 

F e 


Abdo et al. (2011) found that the emission of the main flare 
consisted of five components (see Figure 2 in that work): a nearly 
steady contribution, like a smooth continuation of the plateau 
emission, and four sub-flares of similar duration and energetics. 
Although the uniqueness of such a fit is not statistically assessed, 
it looks like a good empirical description of the data. In the 
framework of the JRGI scenario, such a description is very 
natural. The steady component would be attributed to light 
clouds, continuously ejected by the star, and the four sub- 
flares would correspond to much heavier blobs formed out of 
the blown-up stellar envelop during the initial stage. On the 
other hand, the decomposition of the main flare in four sub- 
flares implies a strict limitation on the variability timescale. The 
flare rise/decay timescales should be longer than the blob light 
crossing time corrected for Doppler boosting. Since the shortest 
variability scale was ~5hr/(l + z rs ) ~ 10 4 s, the maximum 
possible size of the emitting blobs can be estimated as: 

r b « 10 16 r 0 ,i. 5 cm. (18) 

If the jet is magnetically driven, this size constraint can be ex- 
pressed through the mass of the central BH (see Equation (A2)): 

-<0.5 M bi } 9 , (19) 

co 

which is restrictive only in the case of Mbh ,9 1. For 
Mbh ,9 (S 1, the blobs can cover the entire jet without violating 
the causality constraint. 

In summary, if the flare detected with Fermi /LAT was 
produced by an RG entering into the jet, the jet properties should 
satisfy Equations (14), (15), and (17), which correspond to the 
restrictions imposed by the flux level, total energy release, and 
the duration of the plateau stage, respectively. Interestingly, this 
set of equations allows the derivation of a unique solution, which 
can constrain all the key parameters through the value of the 
/7-parameter: 

F 1 - 5 D 1 - 5 

Po = 3 x 10 6 ; 2 7 fc x , 5 erg cm -3 , (20) 

? Z 0,17 


which exceeds the Eddington limit for the mass of the central 
BH M B h ~ 5 x 10 8 M 0 (see Bonnoli et al. 2011 and references 
therein). To assess the feasibility of such a strongly super- 
Eddington jet is beyond the scope of this paper, although 
we note that Lopez-Corredoira & Perucho (2012) presented 
observational evidence indicating that such jets may not be 
uncommon. We also note that the requirement of a high 
luminosity is less severe in the case of small D values and 
that in this limit one can derive a constraint on the central BH 
mass by combining Equations (8) and (23): 

M BH > 7 x 10 8 M q z 3 7 . (26) 

The coherent picture emerging from the jet properties derived 
above suggests that the JRGI scenario can be responsible for 
the flare detected from 3C454.3 for a solution of the problem 
with a reasonable set of model parameters. This solution is 
designed to satisfy the requirements for (1) the total energy; (2) 
the peak luminosity; and (3) the duration of the plateau phase. 
Therefore, some additional observational tests are required to 
prove the feasibility of the suggested scenario. To address 
this issue, we discuss in Section 4 the feasibility of different 
radiation mechanisms to explain the observations under the 
inferred emitter conditions. Also, in Section 3.2, we explore 
whether the stellar wind induced by the RG-jet interaction can 
be responsible for the flux level detected during the plateau 
stage. Recall that, as already shown, the required mass-loss rate 
is well within the limitations imposed by Equation (3). 

Finally, the flare raise time, which is related to the blob 
acceleration timescale (see Equation (5)), can be calculated 
for the obtained solution. Interestingly, in the limit of small D 
values, this timescale appears to be independent of D, the only 
remaining free parameter, and matches closely the detected raise 
time of r r ~ 4.5 hr: 

W(2Ib) -5 hr. (27) 

This agreement can be treated as a cross-check that shows the 
feasibility of the proposed scenario. 


and 


M c = 4M b = 5x10 


f 1-5 ni5 

30 e.max 


r n = 


Fe 2 ' 5 ^ 1 ' 5 


0,17 


F e Z0,17 
Fp max D 


0.5 


S b = 8x 10 


7 0.5 P'1.5 
30 < -0,17 f e 

Z71.5 n0.5 

e,max 


( 21 ) 

( 22 ) 

(23) 


The value of the ram pressure determined by Equation (20) can 
be compared to the value achievable in a magnetically driven 
jet: 

P mdj = 2x 10 4 Zq i 7 ~r^~ erg cm -3 , (24) 

^Edd 


3.2. The Stellar Wind as the Origin of the Plateau Emission 

When the star penetrates the jet, strong perturbations are gen- 
erated in both the jet structure and the star’s external layers. 
Strong shocks and other processes of energy dissipation take 
place in the impacted jet and stellar material, generating con- 
ditions under which nonthermal particles can be accelerated. 
Here, we will not specify the acceleration mechanism, which, 
depending on the acceleration region, may be Fermi I, stochastic 
acceleration, or magnetic reconnection. The accelerated parti- 
cles can either be advected away, while radiating, from the 
perturbed jet region, or penetrate into the stellar atmosphere and 
thermalize heating the ambient plasma. The heating can induce a 
strong wind that will significantly enhance the stellar mass-loss 
rate. A similarly enhanced mass loss can occur as well through 
Compton-heated winds in high-mass X-ray binaries and AGNs 
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due to heating by X-rays from an accretion disk (e.g., Basko 
& Sunyaev 1973; Dorodnitsyn et al. 2008a, 2008b), although 
for the mass-loss rate required to explain the plateau phase (see 
Equation (13)), the accretion disk X-ray flux seems to be too 
small. On the other hand, the nonthermal particles produced in 
the star-jet interaction region can be enough to heat the stellar 
atmosphere. Since even nonthermal particles with low energy 
can effectively heat the atmosphere, the nonthermal emission 
may not be detectable during the heating process. The heating 
energy flux can be estimated as 


F nt = XcPo = 3 x 10 16 / Pq, 6 erg cm 2 s ', (28) 


where Pq, 6 = Po/00 fi ergcm ' 3 ) is the normalized jet ram 
pressure and / = X. d Yd,[ is the nonthermal particle-to-jet power 
ratio within the stellar atmosphere: a combination of the fraction 
of the jet power channeled into nonthermal particles X d and 
the fraction of the nonthermal particles diffusing into the RG 
atmosphere Tdif- In the case of a jet interacting with a heavy 
obstacle, the value of X a could potentially be close to 1. The 
entrainment of particles into the RG atmosphere could also be 
high due to a strongly asymmetric diffusion process: whereas 
the strong magnetic field in the jet may prevent particles from 
crossing back, the weaker magnetic field in the other side would 
stimulate particles to diffuse deep into the stellar atmosphere. 
For all this, the effective value of / may be of the order of 1 . 

The mass flux of a stellar wind excited by nonthermal particles 
of flux F nt can be estimated as: 


B = 10 12 * a_ 12 


1/2 


Gn K 

(CM,) 1 / 2 


= 7x 10 3 a_i 2 P^ 2 2 M S! o /2 xG),6gs 1 cm 2 , 


(29) 


where a is related to the properties of the plasma heating and 
cooling. In the case of X-ray heating, one finds that a ~ 0.03/c 
or ci - \2 ~ 1 (i.e., Basko et al. 1977). However, for nonthermal- 
particle heating, the efficiency of transferring jet energy into 
the stellar wind can be higher because these particles can 
penetrate deeper into the stellar wind. Thus, in this case, a_i 2 
can significantly exceed 1 . The induced mass-loss rate therefore 
will be M w = 7T/2,/? 2 , or 


M « 10 24 a_ 12 <XJ /2 X Po , e g s’ 1 . (30) 

To explain the plateau emission by the interaction of this induced 
stellar wind with the jet, the mass flux from Equation (30) should 
be greater than the one from Equation (13). Adopting the typical 
parameters of the 2010 November flare, this requirement can be 
satisfied if the following condition holds (the function F e here 
characterizes the main flare episode): 

a -i2X > 2 (3D 

Formally, for the adopted RG normalization parameters, the 
required heating efficiency should be rather high: / ~ 0.7 
(accounting for F e < 0.3, as found for conical jets — see 
Appendix B for details). However, a significantly lower heating 
efficiency can be sufficient to fulfill Equation (31) under certain 
circumstances; / could be easily reduced i I D <$/ I (F e ~ D in 
this regime), or assuming that nonthermal particles contribute 
strongly to the pressure of the stellar atmosphere (or__i 2 />> 1), or 
for a RG radius exceeding the fiducial value B* = 100 B 0 . For 
instance, to reduce / to 0. 1, one may adopt either R* = 200 R 0 


or a_i 2 ~ 5 and R., — 100 R e or D ^ 0.05. A combination 
of all these factors could further decrease the required heating 
efficiency to lower values. Therefore, an induced stellar wind 
seems to provide a feasible way to generate the plateau emission 
component. 


4. RADIATION MECHANISMS 
4.1. General Comments 

In this section, we discuss the feasibility of different radiation 
mechanisms for the production site conditions inferred in the 
previous section (see Equations (20)-(23)). The efficiency 
of the radiation channels is determined by the presence of 
nonthermal particles with the required energy and targets. 
Several types of target fields are related to efficient high- 
energy processes: matter in the case of the proton-proton 
or Bremsstrahlung mechanisms, a magnetic field in case of 
synchrotron radiation, and photons for inverse Compton and 
photomeson emission. The obtained solution for the properties 
of the jet, plus observational constraints, allow two of these 
targets to be properly characterized. 

Since the mass and size of the blobs are estimated, one can 
derive the matter density in the production regions: 

n b ~ 10 7 cm" 3 . (32) 


This density estimate allows one to discard the proton-proton 
channel, since the expected cooling time, f pp ~ 10 15 /« s ~ 
10 9 s, is too long. 9 This conclusion also excludes the 
Bremsstrahlung channel, which has a similar cooling time. 

The derived ram pressure of the jet allows the derivation of 
an upper limit on the magnetic field strength: 


Bj < 6 x 10 


F 3 / 4 T) 3/4 

3 1 e,max ^ 
^5/4fci/ 2 ,3/4 
5 *-0.17 


(33) 


In fact, from equipartition arguments, the actual strength of 
the magnetic field should be close to that value; i.e., one may 
take x ~ 1 as the fiducial value for the fraction of the jet 
luminosity carried in electromagnetic form. Also, accounting 
for the relation G?max / F^ 4 ^ 0.7, valid in the range of 
feasible values of the /7-parameter for the conical jet model, 
one can derive a lower limit on the jet magnetic field of 


Bj > 3 x 10 3 r 1/2 * 1/2 2ai7 4 G. (34) 

More detailed calculations of this value are presented in 
Figure 2. In this figure, we also show the value of the 
magnetic field in the comoving system (Bp, which was cal- 
culated using the jet Lorenz factor from Equation (22) as a 
function of the value of D. It can be seen that for the /7-range of 
interest, there is a lower limit on the magnetic field strength: 

Bj = «,r b ' > BjTg 1 

>200r 1/2 >f 1/2 z"i / 7 4 G. (35) 

Regarding the density of the target photon field, it cannot 
be constrained in general. However, in the specific case of the 

9 We note, however, that this limitation does not apply in the case of 

“off-axis” AGNs, since then there are no constraints related to the requirement 

of blob acceleration. In particular, proton-proton interactions were shown to 

be a feasible channel for the interpretation of the TeV emission detected from 

M87 (Barkov et al. 2010, 2012b). 
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Figure 2. Numerical calculations of the parameters related to radiation processes for the case of a conic jet vs. the value of D. Top left panel: jet magnetic field; top 
right panel: jet magnetic field in the comoving frame; bottom left panel: proton synchrotron cooling time (red solid line) and blob acceleration time (blue dashed line), 
both in the comoving frame; bottom right panel: the registered luminosity of the target photon field for the EIC scenario. 

(A color version of this figure is available in the online journal.) 


SSC scenario, it is possible to apply certain limitations. Namely, 
the high-energy component detected by Fermi / LAT exceeds by 
two orders of magnitude (/ rs 100) the flux detected at other 
wavelengths. This finding implies that the energy density of the 
target field should exceed the energy density of the magnetic 
field by the same factor: w' vh = fw' B (these are comoving 
reference frame values). The luminosity of such a field can 
be estimated as 

T target ~ 4jrr b 2 cWp h r£ ~ 4jtr^cfw’ B rl 

~ ^rlcfxP 0 Tl ~ F-' r l fxL Y . (36) 

Since /^target < Ly and F e — 0.5 (T b /r ( |) 2 , one obtains 

xf 2 < 0.5£. (37) 

Thus, the SSC scenario can be realized only in a specific region 
of the jet, where the magnetic field is significantly lower than 
the characteristic one (i.e., x < 10 4 ). The blobs themselves 
may have a magnetic field quite different from the field in the 
jet, and therefore may serve as good sites for SSC. However, this 
assumption involves further complexity, e.g., regarding internal 
shock acceleration within the blobs, and is therefore deferred 
to future studies. Also, we note that the radiation efficiency 
of photomeson production is low (e.g., see the discussion in 


Aharonian et al. 2008). If photomeson production were the 
radiation mechanism, it would imply a very low value of § 
and an uncomfortably high power of the jet. 

We focus now on the following radiation mechanisms: elec- 
tron and proton synchrotron and EIC. The analysis of the syn- 
chrotron channel is straightforward, since the density of the 
target is determined. Indeed, to produce a gamma ray of energy 
Ey.GeV — E y /\ GeV, either an electron or proton of energy 

/ \ 3 / 2 

Eel p = 2 ( ^ ) § I/4 *- 1/4 4 8 7<GeV TeV (38) 

is required. We note that the weak D-dependence in these 
equations is neglected for simplicity. 

To obtain such a high-energy particle, a few conditions should 
be satisfied. The most critical of the conditions are the Hillas 
criterion and the efficiency of the acceleration process. The 
Hillas criterion requires that the gyroradius of the highest energy 
particles is smaller than the size of the acceleration region. Since 
the acceleration region size is determined by the blob cross- 
section, one can derive the following requirements: 

f 3/4 *- 3/4 4Y 7 8 ^ 2 cV < 10 8 ( — ) . (39) 

V m e/p / 

This estimate allows one to illustrate that for the derived 
jet parameters, the Hillas criterion appears to constrain the 


7 


The Astrophysical Journal, 774:113 (15pp), 2013 September 10 


Khangulyan et al. 


acceleration of particles in neither the lepton nor the proton 
synchrotron scenarios. 

Regarding the efficiency of the acceleration process, it is 
convenient to present the nonthermal particle acceleration time 
in a form independent of the specific acceleration mechanism: 


ri R. 


gy 


V Ee/p 

eB' ]C ’ 


(40) 


where E e / p is the energy of particle, r] f 1 is a dimensionless 
parameter characterizing the acceleration process, and R gy is the 
particle gyroradius. This timescale should be compared to the 
dominant cooling time (synchrotron) to determine the highest 
achievable energy: 


9 m 


tsy ~ 4 e*E. 


e/P L 


e/pBj 


> u 


par • 


(41) 


Accounting for the energy of the emitted gamma rays, E Y , one 
can obtain an upper limit on the acceleration time through a 
constraint on the value of i] : 


rj < 5 x 10 3 r b ,! 


Ey 

m e/p c 2 


-l 


(42) 


where T h i = r b /10is the blob Lorentz factor. Since the redshift 
corrected energy of the spectral break appears to be close to 2 
GeV, the leptonic scenario can be realized only if T b > I Or]. 
In the case of proton-synchrotron models, this constraint gets 
significantly relaxed, yielding r] < 200r b . However, given the 
relatively small values of the obtained Lorentz factors (see 
Equation (22)), the derived limits on the acceleration efficiency 
are very tight, especially for the leptonic case (r} — » 1). The 
feasibility of such a scenario should be addressed with a detailed 
acceleration model, which remains beyond the scope of the 
present paper. We note that this limitation can be relaxed by 
adopting a scenario with a highly turbulent magnetic field (see 
more details in Kelner et al. 2013). 

Since the proton synchrotron scenario fulfills the basic 
requirements for the production mechanism, it is worthwhile 
to check whether one can reproduce the observed basic spectral 
properties. Fermi / LAT reported a broken power law spectrum 
with a nearly constant break energy. The change of the power 
law index varies between 0.5 and 1. Whereas index changes 
of 0.5, thought to be cooling breaks, are typical in synchrotron 
scenarios, other values for the power law index change may be 
hard to explain through cooling. On the other hand, one can 
attribute this spectral feature to absorption (as suggested for this 
source by Stern & Poutanen 2011). 

The extrapolation of the synchrotron spectrum below the 
Fermi / LAT range requires a harder slope than at GeV energies, 
n v oc v~ LS or harder, to not overpredict the X-ray flux found by 
Swift. This harder slope implies a break in the photon distribution 
around energies of 10-100 MeV. Such a harder spectrum could 
be associated with the uncooled part of the particle distribution. 
The synchrotron cooling time of protons in the comoving frame 
can be expressed as 


t' . < 10 5 


In Figure 2 (bottom left panel), we show the numerical com- 
putation of these timescales for the case of a conical jet. One 
can see that for the typical model parameters, these time scales 
are comparable. This result allows us to make three important 
conclusions about the considered scenario. Namely, (1) the fast 
cooling regime assumption should not be strongly violated in 
proton-synchrotron models for GeV-emitting protons; (2) one 
should expect a cooling break in the high-energy domain; (3) 
the proton synchrotron mechanism appears to be an efficient 
radiation channel in the Fermi energy range. At lower energies, 
protons could remain uncooled, which may lead to requiring a 
hard particle distribution given the tight energetics of the flares 
in 3C 454.3 plus additional observational constraints (see Sec- 
tion 4.2). It is often argued, based on the comparison of the 
proton cooling and jet dynamical time scales (see, e.g., Sikora 
2011), that proton synchrotron is characterized by a very low 
efficiency as the gamma-ray production mechanism in AGN 
jets. However, as can be seen from Figure 2 (bottom left panel), 
proton synchrotron can produce a high radiation efficiency for 
the inferred jet properties. 

The EIC channel is not strongly constrained in the JRGI 
scenario since the properties of the photon target are not well 
known and can vary within a broad range. Similar to the SSC 
case, the density and luminosity of the target photons can be 
estimated based on the ratio of the synchrotron and inverse 
Compton luminosities. This estimation yields the following 
luminosity: 


Tph.Eic ~ 47rz 2 cwph 

>3 x 10 47 /^ _1 z^5f ergs - '®, (45) 

The value of this luminosity as a function of D is shown in 
Figure 2 (bottom right panel). It is seen that the photon field 
luminosity required for the EIC scenario to work appears to be 
very high (accounting for / ~ 100 10 ), exceeding the reported 
luminosity of the broad-line region (BLR) unless x is very small 
or the interaction region is located far enough from the central 
BH(z » 10 17 cm). 

There is another constraint related to the process of gamma- 
ray absorption through pair creation in the local radiation fields. 
The peaking energy of the radiation component detected with 
Fermi /LAT should be close to 100 MeV. Otherwise, the lower 
energy spectrum should be very hard not to violate the X-ray flux 
detected with Swift (Abdo et al. 201 1). In the emitter comoving 
system, this would correspond to 

E' ph = 10r b - ] ,MeV, (46) 

i.e., relatively close to the electron-positron creation threshold. 
Assuming a photon spectrum oc E ~ 15 below the peak (it may 
be even harder; see Section 4.2), it is possible to compute the 
luminosity transferred to the secondary pairs: 

L sec ~ 5 x 10 48 L 2 50 r K ‘}r- 1 1 5 erg s" 1 , (47) 

where r bi i 5 = r b /10 15 cm is the normalized blob radius. 
Substituting here the solutions obtained for the jet parameters 
(Equations (20)— (23)) in the limit of small D, one finds: 


This timescale should be compared to the timescale character- 
izing the comoving blob acceleration time: 

4c~10 5 4tV- (44) 


L 


sec 


10 


49 - 2.25 

Z 0,11 


ergs 


(48) 


10 A similar consideration to the one done in the context of SSC applies here 
regarding EIC within the cloud, i.e., where potentially x <^\. 
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This energy should be emitted either via synchrotron or inverse 
Compton radiation. The obtained estimate of the proper mag- 
netic field from Equation (35) allows one to estimate the peaking 
energy for the synchrotron channel to be 

7i« se c, sy ~ 10- 2 r 1/2 * 1/2 Zojf eV. (49) 

This component can be constrained by the infrared flux detected 
during the flare (Jorstad et al. 2012). If the dominant channel is 
inverse Compton scattering, the peak energy can be estimated as: 

hco sec , ic ~ 4 (^!L) MeV, (50) 

where 6 eic is the target photon energy. Although this component 
would peak above the X-ray energy band, the lower energy tail 
may provide an important contribution to the reported flux level. 
Also, it is important to note that the inverse Compton photons 
may serve as targets for the absorption of primary gamma rays 
leading to a nonlinear regime of the emission formation. This 
effect, in particular, may make a self-consistent treatment of the 
EIC scenario very complex, and thus we leave this possibility 
for future dedicated studies. Nevertheless, we will account for 
internal gamma-ray absorption in the blob’s own synchrotron 
field and we will comment on radiation reprocessing outside the 
blob in Section 4.2. 

Despite being not very detailed, and order of magnitude, the 
above analysis shows that explaining the 2010 November flare 
detected from 3C454.3 in the framework of the JRGI scenario 
requires a jet with a very high ram pressure. Consequently, a 
magnetic field with a strength not far from equipartition, say a 
factor of 10 below, appears to be too strong for radiation pro- 
cesses other than synchrotron. On the other hand, leptonic syn- 
chrotron emission requires an extremely efficient acceleration 
process with an acceleration parameter t] < 2, which cannot be 
discarded but is perhaps unrealistic. Thus, proton synchrotron 
emission seems the most comfortable radiation channel given 
the restrictions imposed by the data and the JRGI scenario. 
In what follows, we test this process using a simple radiation 
model. 


4.2. Modeling the Lightcurve and the Spectrum 

To check whether the JRGI scenario plus synchrotron radi- 
ation can explain the observations in the case of magnetically 
dominated jets (i.e., x = 1), we have computed the lightcurve of 
the 2010 November flare and the SED for one of its sub-flares. 
The radiation output was assumed to be dominated by proton 
synchrotron; EIC and SSC were neglected due to the strong 
magnetic field (see Section 4.1). 

To derive the lightcurve. Equation (11) was used (see the 
determination of F e in Appendix B). In Figure 3, a computed 
lightcurve that approximately mimics the 2010 November flare 
is presented. The lightcurve has been obtained assuming four 
sub-flares of total (apparent) energy 10 55 erg each, plus a plateau 
component with a luminosity of 2 x 10 49 erg s -1 . For each sub- 
flare, we have adopted D = 0.1. The normalization of the 
lightcurve was determined adopting the following values: the 
Loren tz factor To = 28, the ram pressure Pj = 3x I (/’erg cm 3 , 
blob radius = 2.7 x 10 15 cm, and / = 0.3. These parameters 
imply a minimum jet luminosity of Lj = 2.3 x 10 48 erg s _1 . The 
remaining parameters for the emitter are z o = 1.3 x 10 17 cm, 
and Mb = 1.3 x 10 3 ° g. The corresponding mass of the matter 
lost by the RG to explain the four sub-flares is 5 x 10 3 ° g, not 
far from the upper limit given in Equation (2). 



£ = 0.3. We show four sub-flares (dashed lines), the plateau background (dot- 
dashed line), and the sum of all of them (solid line). The observational data 
points and error bars are obtained from the Fermi /LAT 3 hi' binned count rates 
and photon index using a luminosity distance of Z)l = 5.5 Gpc and assuming a 
pure power law spectrum between 0. 1 and 5 GeV. 

(A color version of this figure is available in the online journal.) 


To calculate the SED, we adopted a spectrum for the injected 
protons Q oc E~ p exp(— E/E cm ) and an homogeneous (one- 
zone) emitter moving toward the observer with Lorentz factor 
Tb = 12. The minimum proton Lorentz factor was taken to 
be equal to the shock Lorentz factor in the blob frame, i.e., 
E mm = r 0 /r b 77ipC 2 . The cutoff energy, E caU was obtained by 
fixing ;/ = 4 x 10 3 (see Equation (40) for details), i.e., a modest 
acceleration efficiency. For the maximum proton energy, i.e., 
how far beyond the cutoff the proton energy is considered, we 
adopted two values: E max = oo and E max = 3 E cut . Regarding 
the latter case, we note that assuming a sharp high-energy 
cut is very natural. The injection spectrum was selected to be 
hard, p = 1, to optimize the required energetics. Interestingly, 
magnetic reconnection, in particular in magnetized jet-cloud 
interactions, has been postulated to provide hard particle spectra 
(see, e.g., Bosch-Ramon 2012 and references therein). The main 
motivation for a hard proton distribution is to avoid the need of 
a VHE budget, although the obtained synchrotron spectrum is 
also consistent with the observed IR and X-ray fluxes, as noted 
in the next paragraph. Significantly steeper particle distributions 
would enhance the required jet power by a factor of a few, but 
could also violate the observational constraints. 

In Figure 4, the SED of a sub-flare is shown. The impact of 
the internal absorption on the gamma-ray spectrum is negligible, 
although the emission of the secondary pairs appear in the energy 
band constrained by optical measurements (Jorstad et al. 2012). 
For the chosen model parameters, the synchrotron secondary 
component goes right through the optical observational con- 
straints, and for slightly higher zo values, the secondary emis- 
sion will be well below the optical points (see Equation (48)). 
Also, we note that the obtained spectrum does not violate the 
X-ray upper limits obtained by Swift. 

To illustrate the impact of external yy -absorption, we have 
introduced a photon field peaking at 40 eV with a luminosity 
of 4 x 10 46 erg s 1 , produced in a ring of radius 10 18 cm 
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Figure 4. Computed SED of the synchrotron emission for a sub-flare of 2010 
November. The thick dashed line shows the intrinsic gamma-ray emission for 
the case of E max = oo. Dotted and dot-dot-dashed lines show the gamma-ray 
spectra corrected for internal absorption only for E max = oo and E max = 3/scuu 
respectively. The thin solid and the dot-dashed lines correspond to the cases 
when absorption is dominated by a blackbody and a monoenergetic photon field, 
respectively. The computed synchrotron SED of the secondary pairs produced 
via internal pair creation is also shown (dotted line). The parameters of the flare 
are the same as in Figure 3. The observational data shown are from Fermi/ LAT, 
Swift (Abdo et al. 2011), and the flux in the R band (Jorstad et al. 2012). 

(A color version of this figure is available in the online journal.) 


at z = 0 around the jet base. Two photon fields have been 
adopted, a blackbody and one represented by a ^-function, to 
simulate the impact of a dominant spectroscopic line. As seen in 
Figure 4, the impact of such an external field can be important. 
The treatment of the secondary emission of the produced pairs 
is beyond the scope of this work. However, we note that the 
ambient magnetic field energy density could be well below that 
of the radiation, allowing electromagnetic cascades to develop, 
effectively increasing the gamma-ray transparency. In fact, if it 
were not the case, the secondary synchrotron emission would 
likely violate the constraints from the Swift data (see, e.g., 
Aharonian et al. 2008; Sitarek & Bednarek 2010; Zacharopoulou 
et al. 2011). Extragalactic background light absorption is not 
relevant in the energy range of interest. 

As seen in Figure 4, the spectral breaks induced by external 
absorption are very sharp and do not fit the Fermi / LAT data 
points. However, the properties of the external photon field can 
vary in a quite broad range, and it is expected that for some 
feasible photon field one could achieve a very good agreement 
between the emission spectrum and Fermi /LAT measurements 
(see, e.g.. Stern & Poutanen 2011). Instead of searching for such 
a field configuration, here we consider an alternative possibility. 
As mentioned above, we assume that the injection spectrum has 
not only an exponential cutoff at E c ut , but also a sharp upper 
limit at E max = 3E cat . As seen in Figure 4, this assumption, 
naturally accounting for the fact that particles cannot have 
arbitrarily high energies in the source, allows one to qualitatively 
model the Fermi /LAT observations without invoking additional 
assumptions regarding external absorption. 

In addition to optical photons, radio emission was also 
detected at the flare epoch and thought to be linked to the 
gamma-ray activity (Jorstad et al. 2012). This radiation is 
strongly sensitive to the details of the flow dynamics, and at this 
stage we will not try to interpret the radio observations. However, 
we note that the energetics involved in gamma-ray production 
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are very large, and the JRGI scenario includes complex MHD 
and radiative processes, so it could easily accommodate the 
presence of a population of radio-emitting electrons at z ^ Zfiare- 

Swift X-rays could be also linked to the JRGI activity. As 
stated in Section 4.1, X-rays may come from secondary pairs 
produced via pair creation, or from a primary population of 
electrons(/positrons). However, as with the radio data, given 
the complexity of the problem we have not tried at this stage 
to explain the X-ray emission contemporaneous with the GeV 
flare. 

5. DISCUSSION 

Observations in the VHE domain put severe constraints on the 
size of the gamma-ray production site, supporting the picture of 
a very compact emitter. The concept of a blob as a compact 
emitter inside jets is often used to interpret the nonthermal 
emission of AGNs. This paradigm is supported by a broad range 
of observations, in particular bright knots moving downstream 
the jets. The origin of such blobs is not clear. It was suggested 
that the blobs form when an external object (e.g., a star or a 
cloud; see Araudo et al. 2010; Barkov et al. 2010 and references 
therein) penetrates the jet and loses its atmosphere due to the 
jet impact. Because of the much larger obstacle density and the 
velocity difference, the intruding matter cannot immediately 
dilute within the jet but has to be first heated and accelerated, 
which in supersonic jets will lead to strong shocks that can 
accelerate particles. Therefore, this scenario naturally presents 
specific characteristics that are hard to introduce in pure jet 
models. For instance, the large blob density makes possible 
the formation of a radiation component typically suppressed 
in AGNs, i.e., emission generated at proton-proton collisions 
(suggested to be behind the TeV emission detected from M87 
by HESS; Barkov et al. 2012b). Moreover, since the size of 
the blobs produced in this scenario is not constrained by the 
gravitational radius of the central BH, it was proposed that 
this process could be responsible for the very rapidly variable 
TeV emission detected front PKS 2155—304 (Barkov et al. 
2012a; Aharonian et al. 2007). The JRGI scenario may have 
quite general applications, but the properties of the emission 
expected within its framework are quite strictly constrained. In 
particular, the available energy budget and expected lightcurve 
shape can be determined with a simple dynamical model, which 
describes the evolution of the blob in the jet. Depending on 
the jet properties and other factors like the viewing angle, the 
dominant emission channel can vary. However, the majority 
of the physical conclusions that can be inferred in the JRGI 
scenario are quite general because they concern the dynamics 
of compact blobs in relativistic jets, independent of the origin 
of the former and the nature of the radiation channel. 

The observations of 3C454.3 with Fermi revealed several 
quite puzzling features, in particular the peculiar lightcurve, 
with a nearly steady plateau phase that was interrupted by 
an exceptionally bright flare. The detected flux corresponds 
to an apparent luminosity of 2 x IIP 0 ergs -1 , which almost 
unavoidably implies a presence of a very powerful jet (see, e.g., 
Bonnoli et al. 2011). In the case of powerful jets, the JRGI 
scenario should proceed in a quite specific way as compared to 
other cases already considered in the literature (Barkov et al. 
2012a). In particular, the mass of the material initially removed 
from the star might be very large, resulting in rather long cloud 
expansion and acceleration timescales, the main flare being 
significantly delayed with respect to the moment of the star’s 
entrance into the jet. The plateau emission would otherwise start 
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just after the jet penetration, and come from the jet crushing the 
lighter clouds ejected from the stellar surface while the star 
travels through the jet. The duration of the plateau phase would 
be determined by the time required by the main cloud to expand 
and accelerate. 

We have studied the lightcurve obtained with Fermi in the 
context of the JRGI scenario, aiming to reproduce three main 
properties of the flare: total energy, maximum luminosity and 
duration of the plateau stage. It was shown that the key properties 
of the jet, i.e., the jet ram pressure (linked to its luminosity) and 
Lorentz factor, as well as the cloud/blob characteristics, i.e., 
mass and cross-section, can be reconstructed as functions of 
the dimensionless parameter D. It was also shown that in the 
limit of small D values, the parameter space is less demanding 
concerning the jet luminosity, and the key characteristics of 
the model saturate at values independent of D , which allows 
conclusive cross-checks of the scenario. In particular, the flare 
rise time appeared to be an independent parameter, with its 
value of 5 hr closely matching the rising time of 4.5 hr obtained 
observationally. Furthermore, it was shown (see Section 3.2) 
that for the inferred jet properties, the jet-induced stellar wind 
can provide a mass-loss rate large enough to generate a steady 
emission component with a luminosity comparable to that of 
the plateau. 

Although the analysis of different radiation channels involves 
additional assumptions regarding the spectrum of the nonther- 
mal particles and the density of the target fields, it was possible 
to show that for magnetic fields not far below equipartition (as 
expected in a magnetically launched jet; see Appendix A for de- 
tails), all the conventional radiation channels can be discarded 
and the emission detected with Fermi can be produced through 
proton synchrotron emission (unless ij — > 1, making electron 
synchrotron also feasible). We note that in this case, the 
emission from pairs created within the blob may also explain 
the reported optical enhancement at the flare epoch. 

If the studied scenario is behind the 2010 November flare from 
3C454.3, it has some important implications on the properties of 
AGN jets. On one hand, the obtained solution implies a super- 
Eddington jet. Although this requirement may appear somewhat 
extreme, given the exceptional properties of the source it cannot 
be ruled out (and in fact super-Eddington jets might be a rather 
common phenomenon; Lopez-Corredoira & Perucho 2012). In 
addition, the obtained solution (Equations (20)-(23)) implies a 
lower limit on the central BH mass that, for a location of the 
interaction region at zo ~ 10 17 cm, is roughly consistent with 
the measured values. Also, adopting a magnetically launched 
jet model, it is possible to infer the mass of the central BH 
as Mbh ,9 = 0.5z 0 17 Tj j 5 , which is again consistent with the 
measurements. 

At the initial stage of the cloud acceleration, the intensity 
of the jet-cloud interaction is the highest, but the associated 
emission is not Doppler boosted and therefore very hard to 
detect in the high-energy regime. Nevertheless, there are other 
manifestations of the IRGI scenario that may be observed. 
In particular, as reported by Leon-Tavares et al. (2013), a 
significant enhancement of the line emission from the source 
was recorded during the flare epoch. The detected flux was 
~2 x 10 45 ergs -1 , implying an ionizing luminosity at a level 
of 10 46 ergs _1 . According to our estimates, the initial blob 
acceleration period should be characterized by luminosities of 
~2 x 10 48 ergs _1 , so only 1% of this luminosity would be 
required to explain the ionization line component. With the 
acceleration of the blob, its radiation gets beamed along the jet, 


thus reducing its impact on the BLR and thereby weakening any 
related line emission. 

Since the duration of the expansion phase determines the 
delay between the onset of the plateau phase and the flare 
itself, it is important to check whether the suggested scenario 
is consistent with other flares registered with Fermi from the 
source, e.g., in 2009 December and 2010 April (Striani et al. 
2010; Ackermann et al. 2010). This issue can be addressed 
through a simple scaling that relates the duration of the plateau 
phase to the total energy released during the active phase: 
r p i oc £ tot (assuming a steady jet, one can obtain from 
Equations (C3) and (C7) that l\) x Tj = const, which yields in 
the scaling dependence). Therefore, for the previous events, with 
energy releases 1-2 orders of magnitude smaller than that of the 
2010 November flare, a rough estimation of the plateau duration 
gives plateau durations between 1.3 and 4 days, consistent with 
observations. 

The rate of the RG penetrating the jet deserves some discus- 
sion. This rate is determined by the stellar density in the close 
vicinity of the central BH, say at distances of ~1 pc. Conven- 
tionally, one adopts a power law dependence of the RG density 
on distance of the sort: p* = po (z / 1 pc ) " . Given that the stellar 
velocity can be estimated as T* ~ hOOMgygr^ 1 ^ 2 kms -1 , one 
can estimate the stellar penetration rate for the distance interval 
(z, z + dz): 

dM * ~ p* V/zT/ 1 dz, (51) 

where we have used z/Tj as the jet radius. 

In the scenario proposed here, the RG penetration takes place 
at a typical height of 10 17 cm. We can therefore roughly derive 
the rate of these events by integrating Equation (51) over the z 
interval ~ (0.3-3) x 10 17 cm: 


M* 


6 x 10 


,-4 


1/2 0-1 

PoA^bh,8 [ 1/2 -a, , 3 — K 

— — X ' dx (pc yr ), (52) 

J o.oi 


where x = z/1 pc is a dimensionless radius. A rather broad 
range of values of the u-parameter is discussed in the literature 
(Bisnovatyi-Kogan et al. 1982; Murphy et al. 1991; Merritt & 
Szell 2006). In particular, Murphy et al. (1991) obtained « = 7/4 
for small stellar densities (~10 6 M 0 pc~ 3 ) and a = 0.5 for high 
densities (~10 8 M 0 pc -3 ). The value of the parameter po can 
be estimated through the mass of the stellar cusp around the 
central BH: p 0 = (1 — a/3)p*, where p* ~ 10 6 ~ 8 M 0 pc~ 3 is 
the averaged (typical) stellar density for z < 1 pc. 

The integral term in Equation (52) spans approximately two 
orders of magnitude, between 0. 1 and 15, for the given range of 
a. We therefore conclude that based on the results obtained 
by Murphy et al. (1991), one should expect a rather weak 
dependence of Equation (52) on the value of a. As a reference, 
we take the high stellar density case with a = 0.5, obtaining: 


M* ~ 6 x 10 


3 -^BH.8 Me 


yr 


(53) 


Given the typical jet Lorentz factor of ~30 and a ~1% fraction 
of RG stars, the above equation corresponds to ~2 events yr -1 
for Mbh, 8 = 1, which is consistent with observations. 

Our study (and before, Barkov et al. 2012a) shows that the 
penetration of a star into the base of the jet in a powerful blazar 
leads to a very characteristic dynamical evolution. In fact, the 
JRGI scenario has quite-constrained physics with basically only 
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one free parameter: the nonthermal efficiency. The specific 
conditions in the jet/blob interaction region will determine, 
as long as particle acceleration occurs, the dominant radiation 
process, but in any case powerful high-emission seems natural 
in the proposed scenario. 
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APPENDIX A 

THE STRUCTURE OF MAGNETICALLY DRIVEN JETS 

Although the process of jet formation is not fully understood, 
recent hydrodynamical studies by different groups have shown 
that the BZ process may be at work in AGNs. These studies 
suggest that, once formed, jets may be magnetically accelerated. 
In this scenario, the jet base is expected to be magnetically 
dominated at z ^ 1 pc (Komissarov et al. 2007; Barkov & 
Komissarov 2008; Beskin 2010). 

During jet propagation, the magnetic field energy can be 
transformed into bulk kinetic energy, with a simple prescription 
for the bulk Lorentz factor at the linear acceleration stage of the 
form (Beskin & Nokhrina 2006): 



(Al) 


where w and r g are the jet and the BH gravitational radius, 
respectively. The opening angle of the jet is expected to be 
0 — w/z ~ 1/Fj (Komissarov et al. 2009). This leads to few 
useful relations between different jet parameters: 


r 2 ^ — 

J 4r„ ’ 


e 2 v ^ ^ — 

z r 2 


4r s . 


(A2) 


It is also useful to express the above relation in the units used 
all through the paper: 

r 2 15 = 0.17(= 0.4 2 )z 0 ,i7M^ 9 , (A3) 


&l\Zo,n — 0.6M B h,9- (A4) 

The jet magnetic field in the comoving frame can also be derived: 


APPENDIX B 

ON BLOB LUMINOSITY IN THE CASE 
OF HEAVY BLOBS 

If the mass of the blob is high, implying a low value of the 
/J-parameter, then blobs can travel a significant distance along 
the jet. Therefore, the jet properties can change, and this should 
have an impact on the emission associated with the blobs. These 
effects can be estimated via the approach suggested by Barkov 
et al. (2012a). The dynamics of the blob are characterized by: 


dt 




p ^ r l 

4 cM h ’ 


(Bl) 


where Pj is the ram pressure of the jet. Since the emission of the 
blob is important only when it moves relativistically, it is safe to 
take z = Zo + ct and the initial condition r h | ; =o = L The above 
equation can be modified to a dimensionless form: 

dg / (ryr 0 ) 2 g 2 \ p 
dy V g 2 (Tj/To) 2 ) (ryro^Po/Pj)’ 


where z = Zoy', wo, To and Po are the jet radius, ram pressure 
and Lorentz factor at z = Zo', respectively, g = Tb/ro; and the 
dimensionless parameter D is determined as 


Pqnrlzq 

4c 2 M h T 3 0 


(B3) 


This approach is nearly identical to the one developed earlier, 
but here one accounts for the possible change of the properties 
of the jet (e.g., for the increase in the jet Lorentz factor). 

If D 1, the blob gets rapidly involved into the jet motion, 
so the solution is identical to the one obtained in Barkov et al. 
(2012a) even if the properties of the jet are changing with z. 
However, if D < 1, the structure of the jet may have some 
influence on the properties of the emission. In particular, for 
purely magnetically driven jets, i.e., with a parabolic shape, the 
above equation can be simplified, since the following relations 
hold: 




co = wo y 1/2 , Pj = Poy \ 


Ij = zo/woy 1 ' 2 = Toy 1 ' 2 . 
Thus, the equation of motion (B2) yields 

dg_ = (y_ _ r;\ D_ 
dy Vg 2 y ) y 2 ' 


(B4) 


with the boundary condition g| v= i = wq/zo- Although this 
equation does not have any analytical solution, an asymptotic 
solution can be derived for D 1 (i.e., the jet structure is 
important): 

ga P = [{wo/zo? + 3D\ny] x /\ (B5) 

For the case of a conical jet (i.e., Tj = const, w = z/T ] and 
Pj = Poy -2 ), the equation of motion reduces to the equation: 


Pc 


2 (L 


1/2 


IZO^lv^ G. 


(A5) 


dg_ 

dy 


= (g~ 


7 D 

g )—, 

y- 


(B6) 
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Figure 5. The derived jet parameters as a function of the Z)-parameter. 
(A color version of this figure is available in the online journal.) 


which allows an analytical solution (for details, see Barkov et al. 
2012a). However, since the analytical solution is rather bulky, 
we use the asymptotic solution: 

gac = [(«oAo) 3 + 3 D(y - 1) /y] 1/3 , (B7) 

which is valid for D <$C 1. Since the models of parabolic and 
conical jets correspond to the most feasible jet model, we will 
consider the realization of the JRGI scenario for these two cases. 
The obtained asymptotic solutions allow one to clearly estimate 
the differences that depend on the specific jet configuration. 

The intensity of the nonthermal emission of the blob, cor- 
rected for Doppler boosting, has the following form: 

L y = ^lJ— ) Tq F e (y) , (B8) 

where F e (y ) = g 4 /(P 0 /Pj)(l/g 2 ~ g 2 /( r j/ r o) 4 ) is the correc- 
tion function. 

The maximum value of the correction function, P cjmix , 
depends on the value of the /7-parameter and on the structure of 
the jet. If D ^ 1 , it saturates at a constant value of /’j. max ^ 0.4 
independent of the jet properties. In the case of D 1, the 
asymptotic solutions show that the maximum of the correction 
function is reached at y ln ~ 1.95, with F eimx ~ 0.8 2/7 2 5 
for a parabolic jet, and for a conical jet, y m ~ 1.33 and 
Te.max ~ 0.46 // 2/! . In the intermediate range, 0. 1 < D < 1, 
Te.max is ~0.5Z>° 5 and ~0.3Z7° 5 for the parabolic and conical 


jets, respectively. The location of the maximum v m has a weak 
dependence on D in this regime; i.e., depending on the structure 
of the jet, the maximum intensity, Fe.max, can vary by a factor 
of ~2 in the case of D < 1 . 

The total energy emitted by a blob can be obtained via the 
integration of the luminosity over the observer time (dr = 
z 0 dy/(2cY 2 )): 

\ 

dyF e j 2 M b/c c 2 T] = 8£ F e M h/c c 2 Y f. (B9) 

This equation defines the function F e . 

For large values of the Z)-parameter, the integration term in 
the above equation can be obtained analytically if one accounts 
for the relation F e = (dg /dy)g 4 / 1): 

00 gmax 

F e = J dyFe^ = J dgg 2 = ^ ax . (B10) 

1 o 

Here we note that, in the case of an accelerating jet, the integral 
determining the total emitted energy is divergent. Formally, this 
is so because g max — > oo. Therefore, the integration should be 
artificially truncated at some point. Physically, this point would 
correspond to the moment when the blob gets homogenized in 
the jet. For large values of Z), the blob speed rapidly reaches the 
jet velocity, and one should adopt g max = 1 . For small I) values, 
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we use ff° dy F e (D / g : ) = I), which is a precise identity for the 
case of a conical jet, and corresponds to a truncation at y ~ 3 
for a parabolic jet. Therefore, for the present order-of-magnitude 
estimates, one can estimate the energy emitted by the blob as 


F e = I dyF e ^r - min ( D, \ 


/■ 


(Bll) 


Since the structure of the jet only imposes uncertainties on 
the order of two, for the order-of-magnitude estimates presented 
in this paper we use the solution obtained for the conical jet 
structure. To use the solution obtained for the parabolic structure 
of the jet, one should properly describe the process of the blob 
homogenization in the jet, which is beyond the scope of this 
work. 


APPENDIX C 

A SOLUTION FOR THE JET PARAMETERS 

One can obtain a simplified analytical solution which allows 
the determination of the jet parameters. Ignoring the weak 
dependence on R t in Equation (17), one obtains: 


^c,30 = 


3 x icr 2 


f Dj.maxT 2 j 5 Sb,32 


e,max A j, 1.5 
^pl 


2 x 10 50 erg 


1 +Zr 

13 days J \ 1.859 

From Equation (15), one can infer the mass of the cloud: 

0.1 ( E y \ 

Mc ’ 3 ° ~~ ? L5 \ 2 x 10 55 erg/ ’ 

and Equation (14) can be re-normalized as 

L 


(Cl) 


(C2) 


Pq = 2 x lO^- 1 F-LS^USF 1 


e,max A 0,1.5°b,32 


2 x 10 50 erg 


erg cm 


(C3) 


Finally, the definition of the D-parameter, Equation (6), provides 
the fourth required equation to determine the parameters Pq, St,, 
M c and To. Here, we assume that the flare was the result of the 
superposition of four sub-flares (i.e., M c = 4Mb): 


D = 


6 x 10 D 0i i 7 

e,max ,30^04.5 


2 x 10 50 erg 


(C4) 


Equations (C1)-(C4) have a unique solution, which character- 
izes the parameters Pq, S b , M c and Tq as functions of D: 


P 0 = 3 x 10' 


rU.5 nl.5 

6 e,max 


F £- 5 %Zo 5 n V 2 x 10 50 erg 


-1.5 


tpi 


13 days 


2 x 10 55 erg 


2.5 


erg cm 


-3 


(C5) 


rl.5 nl.5 

M = 4 Mu — 5 x 10 3Q e,max 

c b J F2.5fc 7 l- 5 

r e 5 ^ 0,17 


2 x 10 55 erg 


2.5 


g. 


2 x 10 5 ° erg 


-1.5 


(C6) 


r 0 = 8 


and 


FeZo, 17 

f n 

1 e.max 


0.5 


2 x 10 5 ° erg 


0.5 


2 x 10 55 erg 


-0.5 


(C7) 


S b = 8 x 10 


0.5 pi. 5 
30 *-0,n r e 


e,max 

E v 


Fq m ax E 0 ' 5 V 2 x 10 50 erg/ \ 13 days 

-1.5 


1.5 


'pi 


2 x 10 55 erg 


cm 2 . 


(C8) 


The values P e ,max and F e depend on the value of the D- 
parameter and the structure of the jet. To illustrate the feasibility 
of the derived solution, in Figure 5 we show the jet parameters 
as functions of the D-parameter. 
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